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Abstract-  The  structure  of  myelinated  axons  is  quite  similar  to 
that  of  a  high-loss  coaxial  cable.  An  electromagnetic  analysis  of 
TM  waves  shows  that  distributed  effects  cannot  be  neglected. 
The  attenuation  and  phase  constants  are  obtained  as  a  function 
of  frequency.  Predicted  finite  wave  delay  in  the  internodal 
segment  approaches  measurements. 
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I.  Introduction 

Electric  stimulation  of  nervous  system  can  restore  motor 
functions  [1].  Electrical  properties  of  nerves  are  ususally 
described  by  means  of  a  cable  model  [2],  [10].  This  model  is 
frequently  simplified,  and  each  internodal  segment  is 
modeled  by  a  lumped  impedance.  This  can  be  a  series  resistor 
[3]-[7],  but  sometimes  parallel  capacitive  impedances  are 
added  [27].  A  lumped  resistor  is  unable  to  model  a  finite 
conduction  velocity  in  the  internodal  segment,  whereas  the 
presence  of  capacitors  modelling  the  myelin  sheath  can 
explain  delays  [27].  Despite  this  fact,  resistor  models  are 
more  frequently  used,  and  the  delay  is  usually  assigned  only 
to  nonlinearities  in  Ranvier  nodes  [15]. 

The  existence  of  travelling  waves  in  internodal  segments 
could  linearly  explain,  at  least  partially,  the  measured 
conduction  speed  [14].  The  objective  of  this  paper  is  to 
determine  whether  under  any  circumstances  a  distributed 
circuit  [28]  covering  the  effects  of  both  the  capacitive  and 
resistive  effects  would  be  more  suitable  than  a  lumped  model. 
To  answer  this  question  we  shall  analyze  an  internodal 
segment  through  which  a  transverse  magnetic  (TM)  wave  is 
guided  [20],  [26].  This  approach  links  with  the  spatially 
distributed  description  of  electric  and  magnetic  fields  found 
in  literature  [7],  [11]-[13],  [18],  [19], 

From  an  electromagnetic  point  of  view,  lumped  element 
models  are  obtained  as  a  quasi-static  approximation  [20]. 
Clark  et  al.  [21]  calculate  the  wave  number  A' as 


For  axoplasm  at  1kHz  with  conductivity  5  mhos/m,  the 
resultant  wave  number  approximates  0.198  rad/m.  In  10  cm., 
the  phase  shift  is  less  than  0.02  radian  or  1.46°,  so  in  [21] 
conclude  the  field  is  quasi-static.  This  conclussion  was  later 
adopted  explicitly  [22]  or  implicitly  in  lumped  element 
models  [3]-[7],  [27].  We  shall  show  that  the  waves  are  not 
defined  by  k,  but  rather  by  the  propagation  constant  /;,  and 
that  quasi-static  approximations  don’t  apply  always,  specially 
for  relatively  high  frequencies. 


II.  Methodology 

For  sake  of  simplicity,  we  shall  consider  an  infinitely  long 
segment,  surrounded  by  an  infinite  extracellular  medium. 
Even  though  the  results  will  not  be  the  same  found  in 
physiological  situations,  qualitative  conclusions  might  be 
useful  to  gain  understanding  into  nerve  conduction.  The 
structure  is  quite  similar  to  a  high-loss  coaxial  cable  (fig.  1). 
Medium  1  is  the  axoplasm,  medium  2  is  the  myelin  sheath, 
and  medium  3  is  the  extracellular  fluid. 


Fig.  1 .  Myelinated  nerve  coaxial  model. 

Axon  radium  is  ai  and  fiber  radium  is  a2. 

We  shall  assume  isotropic,  homogeneous  and  linear  media 
[11]-[13],[16]-[18],  although  the  physiologic  medium  in 
which  the  nerve  is  embedded  is  neither  isotropic  [8]-[9]  nor 
homogeneous  [4].  The  myelin  sheath  is  sometimes  considered 
as  a  perfect  insulator  [3],  [7],  [18].  We  don’t  assume  this 
hypothesis  [27].  All  media  will  be  generic  with  permitivity 
permeability  /J.n,  and  conductivity  c„  where  subscript  i  refers 
to  medium  1,  2  or  3.  Let  the  media  be  charged  due  to  present 
ions.  When  these  ions  come  into  the  axon  at  Ranvier  nodes, 
they  diffuse  [23],  and  thus  a  diffusion  current  appears. 
Therefore,  the  current  density  J  at  any  point  designed  by  r  in 
any  media,  when  an  eletric  field  E  dependent  on  frequency  co 
appears  can  be  expressed  as  [24]: 

j(r,ffl)  =  aE(r,ffl)-X^Vp,(r,®) 

i 

where  the  sum  extends  to  all  present  ions  (Na+ ;K* ;Cl "...)  [23] 
and  Di  is  the  diffusion  constant  of  ion  i.  The  ionic  diffusion 
would  exist  even  if  the  ions  were  uncharged  particles  [24]. 

The  magnetic  field  intensity  H  is  solenoidal  and  may, 
therefore,  be  derived  from  the  curl  of  a  suitable  Hertzian 
vector  potential  function  IIe  [25],  [26]: 
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H  =  j(0£eqV  X  Ue 


(3) 


in  which  Zn  (A.,,/')  is  a  generic  Bessel  function  of  order  n.  By 
definition. 


where 


(4) 


ft) 


Substitution  of  (3)  into  Faraday’s  law  and  subsequent 
integration  leads  to 


E  =  k  2n  e  +  V0 


(5) 


where  k  was  defined  in  (1)  and  0  is  an  arbitrary  scalar 
potential  function.  Using  (1)  in  Ampere-Maxwell’s  law, 


V  X  H  =  jo)£eqE  -  X  DiVpi 

i 

Sustituting  (3)  into  (6)  we  find  that 

V  x  V  x  H  =  k- ne  +  V0  -  X— Vp' 


(6) 


T  j(0£e, 


(7) 


h=4k2-K2 


(12) 


In  these  equations,  hn  =  /3„  -  jan  is  the  propagation 
constant,  inherently  different  to  k.  Phase  constant  is  /)„ 
(rad/m)  and  a„  is  the  attenuation  constant  (Np/m). 
We  shall  proof  that  the  thickness  of  the  axon  myelin  sheath  is 
much  less  than  the  wavelength  of  physiological  signals,  so 
that  the  asymmetric  modes  (dependent  on  0,  and  thus  with 
n  0)  are  evanescent  [25]  and  attenuate  quickly  (P„=0). 
Therefore,  the  main  mode  has  n=0  and  the  fields  will 
correspondingly  be  denoted  E0  and  H0.  The  electric  field  can 
be  obtained  using  (8)  and  (10)  in  (5): 


jh„  dy/n 

E*=-iia*-*+J<afik 


— y  d  — 

-2  ^  ■  dr 


1  ^  1  dp, 

Eog  -  — 

1  r-  dp i 

E0 z  =  «oVo  +  j<W 


(13) 

(14) 

(15) 


Since  0  is  yet  an  arbitrary  function,  we  propose  the 
following  gauge: 


The  magnetic  field  is  found  from  (3): 


<t>  =  v  •  ne 


y  7)  ,  p, 
i  J  (»£,  ,, 


(8) 


Hence  an  homogeneous  wave  equation  is  obtained  for  IE. 
Let’s  adopt  a  cylindrical  coordinates  system  axially  centered 
in  the  axon  ( r,0,z ).  The  wave  equation  reduces  to 

v2nt, +  A:2nE  =  0  (9) 


Despite  the  existence  of  a  difussion  current,  the  election 
of  the  gauge  given  by  (8)  makes  this  equation  to  be  the  same 
as  if  the  media  were  uncharged.  Transverse  magnetic  modes 
may  be  derived  from  ITt,  =  Ke z,  thus  assuming  Tle  have  no 
other  components  but  the  axial.  The  solution  to  wave 
equation  (9)  is  [25]: 


n 


n=0  \ 


(10) 


where  an  is  an  amplitude  constant  and  i//n  is  a  is  a  space  and 
frequency  dependent  function  given  by 


Hn. 


H0r  =  0 

_Jr_  dy^ 

poA2  U°  dr 

H0z  =  0 


(16) 

(17) 

(18) 


Note  the  absence  of  axial  component  in  the  magnetic  field 
in  this  propagation  mode,  despite  the  diffusion  current.  This 
is  therefore  a  TM  mode  [26].  At  r=0  the  field  must  be  finite 
so  Z,''1  =  J0,  with  J0  the  Bessel  function  of  the  first  kind  and 
order  0.  In  the  external  fluid,  the  field  can  not  increase  with 
distance  and  should  behave  like  a  travelling  wave; 
therefore  Z^J)  =  H {0,) ,  where  H0  is  the  Hankel  function,  fn  the 
myelin  sheath,  Z(02)  =  J0  +c„Ng,  with  N0  the  Bessel  function 
of  the  second  kind  and  order  0,  and  c0  a  constant  to  be 
determined. 

Boundary  conditions  for  the  continuity  of  E:  and  Hq  at 
r=ai  y  r=a2  lead  to  an  inhomogeneous  equations  set  for  a(01] , 

c/q2>  ,  c0a(02>  and  a(g3)  ■  For  the  system  to  have  more  than  a 
single  solution,  and  allow  waveguiding  happen  independently 
of  charge  gradient,  the  determinant  of  this  set  of  equations 
must  vanish.  The  determinantal  equation  is  the  same  that  can 
be  found  for  uncharged  media  [25]: 


V,  =  Z„(X„r)e-tMe-*- 


(11) 
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AjAf  \k;  ,/0(/l1a1)/V1(A,<5[1) 

k2k\  —  \h1JQ(hia^)Jl(?i1a^  (19) 

_  ^Hl1){^a2)N0(^ai)- Q 

KklH\'}  (  V2  Vo  ( Vl  )  -  \klH«]  (V2  Vl  (  V2  ) 

The  roots  of  equation  (9)  together  with  (12)  are  the 
propagation  constants  h0m,  with  m=l,2,3.... 

III.  Results 

The  electrical  properties  of  media  shown  in  Fig.  1  can  be 
found  in  Table  I.  We  are  interested  only  in  hoi,  since  higher 
modes  present  stronger  attenuation.  Figs.  2  and  3  show  the 
the  phase  constant  /3  and  the  attenuation  constant  a  as  a 
function  of  frequency.  Frequency  values  have  been  chosen  so 
as  to  cover  the  whole  spectrum  of  physiological  signals. 

Signal  velocity  can  be  calculated  as  co/j3  (fig.  4). 
Wavelenghts  are  determined  as  2n/p  (see  Table  II). 

Table i 


ELECTRICAL  PROPERTIES  OF  MYELINATED  NERVE  MEDIA 


Parameter 

Symbol 

Value 

Unit 

Reference 

Axoplasm  Conductivity 

o. 

0.909 

S/m 

[29] 

Myelin  Conductivity 

a2 

3.45e-6 

S/m 

[30] 

External  Conductivity 

a3 

0.333 

S/m 

[31] 

Axoplasm  Permitivity 

El 

£o 

F/m 

Myelin  Permitivity3 

£2 

25’5eo 

F/rn 

External  Permitivity 

£3 

£o 

F/m 

Axon  Radium” 

ai 

6.2 

pm 

Fiber  Radium 

a2 

10 

pm 

[23] 

“Derived  from  coaxial  capacity  5e-3  gF/cnf  [10],  [24] 
““Using  2a,  =  0.8  x  2a2  -  1.8  [32] 
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Fig.  2.  Phase  constant  as  a  function  of  frequency 


Propagation  Mode:  TM 


Frequency  (Hz) 


Fig.  3.  Attenuation  constant  as  a  function  of  frequency. 


Propagation  Mode:  TM 


Frequency  (Hz) 

Fig.  4.  Wave  velocity  as  a  function  of  frequency. 

IV.  Discussion 

Phase  constant  is,  at  1kHz,  higher  than  200rad/m  (125.6° 
in  10cm).  This  value  is  far  away  from  0.198  rad/m  pointed 
out  by  Clark  et  al.  [21],  [22]  and  shows  the  need  for  a 
distributed  parameter  model.  Note  that,  as  one  could  expect, 
this  is  more  true  as  frequency  increases.  Therefore,  the  nerve 
behaviour  is  not  always  quasi-static.  Wavenumber  and 
propagation  constant  are  the  same  when  the  wave  is  TEM 
[26],  but  we  are  dealing  with  TM  modes.  Attenuation 
constant  is  more  or  less  constant  at  frequencies  below  250 
Hz,  and  then  increases  linearly.  For  a  typical  internodal 
distance  of  1  mm  [23],  low  frequency  signals  attenuate  up  to 
0.28  times  their  value,  despite  their  frequency.  At  higher 
frequencies,  attenuation  grows  approximately  at  a  3 
dB/decade  rate. 

The  order  of  magnitude  of  concerning  wavelengths  are  a 
thousand  times  greater  than  the  thickness  of  the  myelin 
sheath,  validating  a  previous  assumption. 
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Signal  velocity  is  independent  of  frequency  up  to  500  Hz. 
Values  are  similar  to  those  measured  in  [33]  (25-41  m/s), 
[14].  Therefore  delay  is  not  only  due  to  Ranvier  nodes,  but 
also  to  internodal  segments.  This  result  is  in  agreement  with 
[27],  but  no  we  have  used  no  lumped  element. 

Equation  (2)  is  an  extension  of  Ohm’s  law.  Lumped 
element  circuits  neglect  the  diffusion  current,  taking  into 
account  only  the  drift  current  due  to  the  electric  field.  But 
ionic  diffusion  would  be  present  even  in  absence  of  axial 
electric  field. 

Table  ii 

_ WAVELENGTHS  IN  MYELINATED  NERVES _ 

Frequency  (Hz)  Wavelength  (mm) 


50 

396 

100 

216 

500 

46 

1000 

24 

2500 

10 

5000 

6 

V.  Conclusion 

Whereas  lumped  circuit  models  are  much  easier  to  deal 
with,  the  quasi-static  hypothesis  doesn’t  hold  for  relatively 
high  frequencies  and  a  distributed  model  might  be  more 
adequate.  TM  waves  linearly  model  attenuation  and  non-zero 
delay  in  the  internodal  segment,  and  could  help  to  gain  more 
insight  into  neural  artificial  stimulation  problems.  Further 
research  is  required  to  analyze  other  propagation  modes. 
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